Coagulation and fragmentation dynamics of inertial particles 
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Inertial particles suspended in many natural and industrial flows undergo coagulation upon col- 
lisions and fragmentation if their size becomes too large or if they experience large shear. Here we 
study this coagulation-fragmentation process in time-periodic incompressible flows. We find that 
this process approaches an asymptotic, dynamical steady state where the average number of parti- 
cles of each size is roughly constant. We compare the steady-state size distributions corresponding 
to two fragmentation mechanisms and for different flows and find that the steady state is mostly 
independent of the coagulation process. While collision rates determine the transient behavior, 
fragmentation determines the steady state. For example, for fragmentation due to shear, flows that 
have very different local particle concentrations can result in similar particle size distributions if the 
temporal or spatial variation of shear forces is similar. 

PACS numbers: 05.45.-a, 47.52. +j, 47.53. +n 



I. INTRODUCTION 

The dynamics of inertial particles in fluid flows plays 
an important role in many natural and industrial con- 
texts and has received an increasing interest in recent 
years. Questions where inertial particles play an impor- 
tant role are ubiquitous in biology, chemistry, oceanogra- 
phy, astro- and geophysics. Recent works from dynami- 
cal systems [1, 2, 3, 4, 5, 6, 7, 8] to atmospheric science 
[9, 10, 11] and turbulence [12, 13, 14] have added greatly 
to the understanding of these phenomena. Almost all 
these works have been devoted to the dynamics of non- 
interacting inertial particles. A major reason for this is 
that this problem is already very rich, displaying features 
yet to be understood in their full complexity, as inho- 
mogeneous spatial distributions [15, 16] and multivalued 
velocity fields [12, 17, 18]. Interestingly, these very same 
features yield an increased rate of collisions [13], the con- 
sequences of which are in most cases not explicitly taken 
into account. Typically one assumes a dilute regime and 
fully neglects the collisions. In some other works, one 
keeps track of the collisions numerically without actually 
addressing the outcome of such events {ghost collisions) 
[13, 19, 20]. To our knowledge, only very recent works 
have addressed effects of collisions on the dynamics of 
inertial particles [21, 22, 23]. In Ref. [22], we have re- 
ported our first results on the dynamics of inertial par- 
ticles coagulating [48] upon collisions and fragmenting 
under certain conditions. In Ref. [23], the authors con- 
sidered elastic collisions in a monodisperse system and 
pointed out the existence of bursts in the spread of the 
particles out of the attractors of the purely advectivc dy- 
namics. In Ref. [21] the authors treated coagulation and 
shear-fragmentation of dust particles in an astrophysical 
context. There small dust particles can grow into larger 
fractal clusters due to turbulent collisions. 

In this paper we extend the work of Ref. [22] to differ- 
ent flows and to a broadened parameter set. Our motiva- 



tion lies primarily on natural phenomena such as the col- 
lisional growth of cloud droplets [24], sediments in lakes 
and rivers, and marine snow in the ocean [25]. Here we 
focus on the description of spherical droplets, i.e. we do 
not take into account any fractal structures that often 
appear in sediments or marine snow. 

Our main result is that coagulation and fragmenta- 
tion dominate the behavior of different time spans of the 
process and fragmentation rather than coagulation is the 
dominating process for the steady state size distribution. 
For different flows, the collision rates between inertial 
particles can be very different, leading to great changes 
in the coagulation of particles. While this might be an 
important effect for transient processes, such as the ini- 
tiation of rain in clouds, it will turn out that for the 
steady state fragmentation plays a much greater role. In 
fact, the steady state size distributions are mainly deter- 
mined by the fragmentation process. The specific flow 
structure is only relevant for the steady state size distri- 
bution when it directly affects the fragmentation process. 
This can for example be the case for fragmentation due 
to shear when the spatial and temporal variation of the 
shear is very different for two flows. 

We consider fragmentation to be of two possible ori- 
gins. First, particles break up if their size exceeds a cer- 
tain maximum allowed size. This is motivated by the hy- 
drodynamical instability of large water drops (e.g. cloud 
drops) settling due to gravity [26]. Second, particles frag- 
ment if the shear forces due to the fluid flow are suffi- 
ciently large. This mechanism has been reported to be 
the dominant one in the case of marine aggregates [27]. 

At a first glance, one might be tempted to pursue a 
field-theoretical approach, in the framework of which one 
treats the problem of particle motion as a multiphase 
flow and then applies the Smoluchowski equation [28] to 
model coagulation and fragmentation for the particle dis- 
tribution. However, the inertial particle dynamics is dis- 
sipative and contracts to an attractor in a 2(i-dimcnsional 
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phase space, where d is the spatial dimension of the flow. 
This attractor can be folded in phase space, meaning that 
the particle velocity may take on several values even at 
the same location. Due to the presence of such 'caustics' 
[12, 18, 29], a field-theoretical approach cannot be well 
founded. Therefore a study based on an individual track- 
ing of the particles, as the one presented here, becomes 
necessary. 

Here we consider the fluid flow to be spatially smooth 
and to have a single macroscopic time scale. We are mo- 
tivated by flows having coherent (e.g. convective) struc- 
tures on length scales much larger than the ones at which 
turbulence plays a major role. The effect of turbulence 
can then be taken into account as a stochastic perturba- 
tion described by an eddy diffusivity [30] at small scales. 
For simplicity, we neglect this small scale noise in the 
present work and focus only on the large scale motion of 
the fluid. 

We study the dynamics of the system formed by the 
transported inertial particles undergoing coagulation and 
fragmentation in three different fluid flows, as described 
in Section II. We find that the system tends to approach 
a steady state where several size classes coexist (Section 
III) . The average number of particles in each size class is 
roughly constant with a mild periodic time dependence 
— with a period identical to the one of the advecting 
fluid flow. The distribution of particles as well as the 
mean average size in the steady state depends on the 
type of fragmentation mechanism taking place. First, 
when fragmentation occurs solely due to particles ex- 
ceeding a maximum allowed size, the distribution is in 
general quite broad. Second, for fragmentation occur- 
ring also under sufficiently large shear, the distributions 
typically decay exponentially beyond a certain size class. 
The distributions depend on the fluid flow for both types 
of fragmentation. However, for shear fragmentation the 
differences are very small as long as the variation of the 
fluid shear over time is qualitatively similar. 

In the case of shear fragmentation, we derive a scaling 
relation for the average size class in the steady state as a 
function of the coagulate strength parameter 7. Finally, 
we show that our results are robust with respect to the 
total mass of particles, the number of allowed size classes 
and the initial particle size distribution. Also, in the case 
of shear fragmentation the size distribution in the steady 
state has a scaled functional form which does not depend 
on the coagulate strength 7. 



II. COAGULATION AND FRAGMENTATION 
MODEL 

A. Dynamics of Inertial Particles 

First, we present the equations of motion for the mo- 
tion of finite size particles that will be used here. For sim- 
plicity we consider heavy spherical aerosols, i.e. particles 
much denser than the ambient fluid and assume that the 



difference between their velocity x and the fluid velocity 
u = u(x(i),t) at the same position is sufficiently small 
so that the drag force is proportional to this difference 
(Stokes drag). The dimensionless form of the governing 
equation for the path x(t) = (xi{t), X2{t)) of the center 
of mass for such heavy aerosols subjected to drag and 
gravity reads in this case as [31, 32, 33]: 

x=-(u(x(0,t)-i-iyn), (1) 

T 

where n is a unit vector pointing upwards in the vertical 
direction. Throughout this paper we consider the vertical 
direction along the axis X2. Under the assumption that 
the density ratio p f / pp <^ 1 , the particle response time t 
can be written in terms of the density pp of the particle, 
the radius r of the aerosols, the fiuids dynamic viscos- 
ity 7], and the characteristic length L and characteristic 
velocity U of the flow as t = {2r'^ ppU)/{9riL). Wc note 
that the response time t is nothing but the Stokes num- 
ber which can be written in our case as r = Tp/T where 
Tp is the particle's dimensional Stokesian relaxation time 
and T is the characteristic time of the flow. The dimen- 
sionless settling velocity in a medium at rest is given by 
W = {2r^ppg)/{9r]U). Note that W/t is independent of 
the particle radius r. 

Every particle produces perturbations in the fiow that 
decay at least inversely proportional to the distance from 
the particle [34, 35]. Here we assume a dilute regime, 
where the local concentration of particles is low enough, 
so that particle-particle interaction can be neglected [36] 
unless particles come into direct contact. 

The assumption that the particle radii a are small also 
means that the feedback from the particle motion on the 
fiow will be small as well [33] and is therefore neglected 
in the following. 

B. Coagulation 

Second, we present a model for the coagulation of finite 
size particles. 

The smallest particles considered will be called pri- 
mary particles. These primary particles can combine 
to form larger particles, called coagulates. Coagulation 
takes place upon collision. All particles are assumed to 
consist of an integer number of these primary particles, 
i.e. the primary particles can never be broken up. The 
number a of primary particles in a coagulate is called the 
size class index. We consider n different size classes, i.e. 
coagulates can consist of a maximum of n primary parti- 
cles. A coagulate of size class a has a radius 7'q ~ a^^^ri, 
where ri is the radius of the primary particles. The re- 
sponse time is Tq, = {ra/ri)^Ti = c?^'^t\ and the settling 
velocity in stiff fluid is Wa = c?I^W\. Here n and W\ 
are the response time and the settling velocity for the 
primary particles, respectively. The largest coagulates 
therefore have a radius r„ = v}l^r\. We note that parti- 
cles of different sizes have different parameters Tq and Wa 
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and therefore follow the flow with different parameters in 
the equation of motion (1). 

We define a collision of two particles if the centers 
of the particles, say of radius and Vj, come closer 
than a distance d — ri + rj. In that case the parti- 
cles coagulate and form a larger particle. Mass con- 
servation requires the radius of the new particle to be 
'''new = + ^'j- For the size class index this implies a 
linear rule, anew = ct? + ctj, which determines the new 

2/3 

response time and settling velocity via t^^^^^ = anewTi 

2/3 

and Wa^^^ = anewWi, respectively. 

The velocity of the new particle follows from momen- 
tum conservation. The position of the new coagulate is 
the center of gravity of the two old particles. 



C. Fragmentation 

Third, we present a model for the fragmentation of 
particles. Primary particles cannot be broken up. In the 
following, we will compare two different fragmentation 
rules. 

(i) Size-limiting fragmentation: If a particle becomes 
larger than the maximum radius a„, it is broken up 
into two smaller fragments (binary fragmentation) whose 
radii are chosen randomly, from a uniform distribution 
between ai and half the original radius. If any fragment 
is larger than a„ this process is repeated, until no frag- 
ment exceeds a„. 

(ii) Shear fragmentation takes place when the hydro- 
dynamical force Fhyd acting on the particle exceeds the 
forces Fcoag holding the coagulate together by a certain 
factor. The criterion for breakup can therefore be ex- 
pressed as 

Fhyd/Fcoag > 7 (2) 

where 7 is a constant. 

The hydrodynamical force in this case is proportional 
to the local velocity gradients in the flow. It is expected 
that larger particles are more likely to break-up, there- 
fore the critical force required for fragmentation should 
decrease with the coagulate size. For liquid spherical 
particles (drops) in the size range where viscous forces 
dominate, Taylor [37] and later Delichatsios [38] derived 
an expression for the critical velocity difference Au across 
the drop required for breakup. Under the condition that 
the characteristic time of drop deformations is small com- 
pared to the time where this velocity gradient occurs, we 
rewrite this condition with a single parameter as 

-=7(^)=7«-/^ (3) 

where 7 is a constant, the coagulate strength parameter 
(the same quantity is called stickiness in [22] ) . The radius 
has been normalized with the radius of a primary parti- 
cle. If the maximum velocity difference across the radius 



of the drop exceeds the threshold value given by Eq. (3), 
the particle is broken up into two smaller fragments (bi- 
nary fragmentation) in the same way as for size-limiting 
fragmentation. 

At the instant of both coagulation and fragmentation 
there is a sudden change in the dynamics: the number of 
particles changes in 2 or 3 among the n available dynam- 
ical systems defined by the size classes. 



D. Fluid flows 

For convenience, we treat the case where the fluid flow 
depends only on two coordinates, i.e. we study a three- 
dimensional flow where the velocity in the third direc- 
tion is negligible compared to the other two velocities. 
This can then be represented by a two-dimensional flow. 
Therefore, the phase space of the particles dynamics is 
four-dimensional. We choose three simple paradigmatic 
flow situations with different characteristics to indicate 
the generality of our results. 

All flow domains are spatially periodic, with a charac- 
teristic length L. More specifically, the flows are (a) a 
convection cell flow with moving vortex centers (in the 
following referred to as the moving convection flow), (b) 
a convection cell flow with fixed vortices (referred to as 
the fixed convection flow), (c) a sinusoidal shear flow. 

The two convection cell flows (a) and (b) consist of 
a regular pattern of vortices, or roll cells. Flow (b) 
was first introduced by Chandrasekhar [39] as a solution 
to the Rayleigh-Bcnard problem and since then it has 
been used in the context of different theoretical stud- 
ies [1, 40, 41, 42]. The moving convection flow (a) is a 
slightly modified version, with moving vortex centers, to 
yield a more realistic chaotic regime for the particle mo- 
tion. Convection flows are chosen because they contain 
vortices (convection cells) and uprising/sinking regions, 
which are characteristic features of realistic flows often 
found in nature. The flows are defined by the velocity 
field as 

(a) moving convection flow 

/ ^-\ n , ; • / +M sin(27rx"i) cos(27rf2) ^ 

u(xi,a;2,t) = [l+fcism(c.ii)](^_eos(27rfi)sin(27r^";) ) ' 

(4) 

where Xi = Xi + k2 sm(uj2t) and X2 = X2 + k2 coa{uj2t). 
The parameters ki = 2.72 and uji = tt are the amplitude 
and the frequency of the periodic forcing of the flow, re- 
spectively. ^2 = l/(27r) and UJ2 ~ 7r/4 determine the 
amplitude and the frequency of the periodic motion of 
the centers of the vortices in the flow. The period of the 
flow is T = 2 and the characteristic length and charac- 
teristic velocity are L = 1 and U ^ 1. 

(b) fixed convection flow with the same equation for the 
flow as in (a), but with ^2 = 0. 

The sinusoidal shear flow (c) consists of alternating 
horizontal and vertical velocity components, where each 
velocity component consists of two plateaus in time. It 
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was introduced in Refs. [43, 44] and has been used many 
times in chaotic advection studies. Here we consider a 
time-continuous version (see [45]) defined by: 
(c) sinusoidal shear flow 

+ tanh(/?sin(2^t))) sin(27ra;2) \ 
ui^a;i,a;2,rj-u.o^^^^_^^^^j^^^_^.^^27rt)))sin(27ra;i) J ' 

(5) 

where the parameter (3 describes how rapidly the transi- 
tion between two values, a zero and a nonzero velocity, 
takes place for each velocity component. The typically 
used value (3 = 20/t: corresponds to a very rapid transi- 
tion. 

The period of the flow is T = 1 and the characteristic 
length and characteristic velocity are L = 1 and U = 1. 
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FIG. 1: Maximum of the velocity gradient grad u{xi,X2,t) 
vs. time for the two different flows described by Eqs. (4) and 
(5). 

The fluid flows are laminar and time-periodic, but the 
dynamics of the inertial particles moving in these flows 
can be chaotic. 

To emphasize the difference between the flows, Fig. 1 
shows the maximum of the velocity gradient vs. time 
for the convection flows (there is no difference between 
the moving and the flxed convection flow) and the sinu- 
soidal shear flow. The difference in magnitude and also in 
the temporal evolution between these two flows is clearly 
visible, indicating a possibly very different behavior with 
respect to shear fragmentation. 



E. Numerical implementation 

After presenting the model, we describe some details 
about the implementation. In the bulk of the paper we 
consider n = 30 size classes. The primary particles con- 
sidered here have dimensionless radius ri = 5/30^^'^ x 
10~^, response time ti = 1/55 and settling velocity 
Wi = 3.2ri. As initial condition we take N{t = 0) = 10^ 
particles of the smallest size and no larger particles. Fur- 
thermore, particles are uniformly distributed over the 
1x1 unit cell of the configuration space. The initial 
particle velocity matches that of the fluid at their posi- 
tion in all cases. 

The simulation is based on the following ingredients: 



a) All particles move in the flow over some time step 
dt according to Eq. (1). This integration time step dt 
needs to be chosen small enough to allow for the detec- 
tion of every collision. After each time step dt there is 
an interaction between particles in the form of coagula- 
tion if they are too close to each other. Our experience 
shows that a choice dt = T/20 is sufficiently small for the 
conditions considered here. 

Because of the spatial periodicity of the fiow, the parti- 
cle dynamics is folded back onto the 1x1 unit cell, using 
periodic boundary conditions (see e.g. [40, 42]). 

b) Particles coagulate if their distance is smaller than 
the sum of their radii. Computationally, the coagulation 
process is the most costly component of the simulation. 
Here a link-cell algorithm [46] is used to compute the 
distance between particles, which scales as 0{N) and is 
thus much faster than simply summing over all particles. 

c) Coagulates can fragment cither due to size-limiting 
fragmentation or due to shear fragmentation. 

1. Size-limiting fragmentation: If the coagulate size 
a exceeds the predefined maximum size, which is 
in the following fixed at n = 30 unless mentioned 
otherwise, the coagulate is broken up. 

2. Shear fragmentation: If the shear at the position of 
the coagulate exceeds a critical value, determined 
by Eq. (3) the coagulate breaks up. Due to the 
symmetry of the flows chosen here, the maximum 
velocity difference is always in the direction of one 
of the coordinate axes, therefore only these values 
have to be calculated. Shear fragmentation is al- 
ways applied together with size-limiting fragmen- 
tation to keep the maximum number of occurring 
size classes fixed at n. 

Whatever rule is applied, the result is the reversed pro- 
cess of coagulation: two new particles are formed from an 
old one with the size class indices: ai,new+ctj,new = ctoid- 
As indicated earlier, at^new can take on any value between 
Qfi and an/2 with equal probability. The centers of the 
new particles are placed along a line segment in a ran- 
dom direction so that the distance d between the particle 
centers equals the sum of their radii, i.e. d = ri+rj, and 
the center of mass remains unchanged. Momentum is 
conserved. For simplicity we assume that the new parti- 
cles have the same velocity as the old one. 

III. SIMULATION RESULTS 

In this section we show simulation results using the 
model described above and compare the influence of the 
different flows, and the effect of size-limiting fragmenta- 
tion and shear fragmentation. 

Before presenting any results for the complete model, 
it is worth showing the attractors for the non-interacting 
problem in the different flows. Fig. 2 shows the strobo- 
scopic section (taken with the period T of the flow) of the 
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attractors for the flows (a)-(c) for one specific size class 
projected onto the plane of the coordinates. The figure 
illustrates the difference in the geometric properties of 
the particles dynamics in the different flows. 



(a) (b) (c) 




FIG. 2: Stroboscopic sections of the attractors of Eq. (1) 
projected onto the configuration space for particles for size 
class a = 16 for: (a) moving convection flow, (b) fixed con- 
vection fiow and (c) sinusoidal shear flow. The difference in 
the degree of clustering in the three fiows is clearly visible. 

For the moving convection flow and the sinusoidal 
shear flow the degree of clustering of the particles in 
the attractors, quantified by their fractal dimension, de- 
creases monotonically with the size class. The parameter 
region is chosen in such a way that the attractors are ei- 
ther area filling or fractal with dimension smaller than 2, 
which we consider to be closer to a realistic situation than 
for example fixed point attractors. For the fixed convec- 
tion flow, the particles tend to cluster on a quasiperiodic 
attractor. This leads to a much larger collision rate than 
in the other two cases. 

In all fiows we find convergence to an asymptotic 
steady state. Initially, coagulation leads to a fast increase 
in the average particle size class, independent of the frag- 
mentation rules. Then fragmentation sets in and a bal- 
ance between coagulation and fragmentation is reached, 
with an asymptotic average coagulate size 

t+T 

aoo = lim — / ds {a{s)) (6) 

t^oo ± J 

t 

that depends on the fragmentation rule and the differ- 
ent flows. The average {a{s)) is taken over the coagulate 
sizes at time s and Eq. (6) corresponds to time-averaging 
(a(s)) over one time period of the fiow to remove the pe- 
riodicity. For the transient behavior of the coagulation- 
fragmentation process the geometric properties, in par- 
ticular the degree of clustering of the particles, are very 
important since they affect the time scales of the tran- 
sients. This is not the case for the steady state. This 
becomes very clear when looking at the size distributions 
of coagulates in steady state for the different flows and 
fragmentation mechanisms (Fig. 3). One might expect 
that the steady state for flow (b) is very different from 
the other two cases, due to the large difference in parti- 
cle clustering. Our results show however, that the steady 



state of the particle dynamics for fiows (a) and (b) are 
almost identical, while fiow (c) produces different results. 
For example for both fiows (a) and (b) the size distribu- 
tion has a long tail towards larger size classes that decays 
exponentially for shear fragmentation. Similar exponen- 
tial tails for the particle size distribution have been found 
in observations of aggregates in the ocean (see e.g. [47]). 

By contrast, for the sinusoidal shear fiow the size distri- 
bution has two peaks and then drops off sharply towards 
zero beyond the second peak. In this case the size dis- 
tribution for shear fragmentation is almost identical to 
that of size-limiting fragmentation, but for a lower value 
of oiraax- This is duc to the fact that in the sinusoidal 
shear flow for the chosen parameters the shear a particle 
experiences is almost constant over time and space, ex- 
cept for a small 'dip' every half period. This very narrow 
distribution of the shear is very similar to having a sin- 
gle maximum stable size, as is the case for size-limiting 
fragmentation. In this case a value of 7 = 17 for shear 
fragmentation corresponds to a value of Umax = 20. 

We have also checked the size distributions in sub- 
regions of the flows. We found, that as long as the 
number of particles in the subregions allowed sufficiently 
good statistics, the normalized size distributions coin- 
cided with the global distributions. This means that 
there is no signiflcant spatial dependence of the size dis- 
tribution. 

While the speciflc shape of the size distributions found 
may not be very general, both because of the limitation to 
only thirty size classes and the very simplifled flows, our 
approach illustrates clearly that the geometric properties 
of the particle motion related to preferential concentra- 
tion are not the most relevant ones for the steady state 
distributions. Instead, the most important effect for the 
steady state of the particles seems to be the fragmenta- 
tion process. 

Size-limiting fragmentation is the same in all flows, as 
it does not depend on properties of the fiow. In this 
case the differences in the size distributions are small 
for the different flows. This indicates that the size dis- 
tributions are mainly determined by the fragmentation 
process, because the flow specific differences, e.g. differ- 
ences in coagulation rates, do not affect the shape of the 
size distribution. However, when fragmentation depends 
on shear, the different flows produce very different size 
distributions. This difference is mainly due to the dif- 
ferent properties of the shear forces in the flow (Fig. 1) 
which lead to differences in fragmentation. To see that 
it does not depend e.g. on the detailed characteristics of 
the particle motion or on the different collision rates, we 
can adjust the flow parameters. Decreasing the value of 
the parameter (3 for the sinusoidal shear flow to a much 
smaller one, e.g. /3 = 7r/20 greatly increases the variation 
in the shear forces over space and time. We obtain two 
sinusoidal peaks per period for the shear forces, similar to 
what happens in the convection fiow (cf. Fig. 1), except 
that for the sinusoidal shear fiow both peaks have the 
same height. For an appropriate choice of 7, so that the 
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FIG. 3: Histogram of the particle size distribution in steady 
state for size-limiting fragmentation and shear fragmentation 
with (a) moving convection cell flow, (with 7 = 50 for shear 
fragmentation), (b) fixed convection cell flow, (with 7 = 45 for 
shear fragmentation), (c) sinusoidal shear flow (with 7 = 17 
for shear fragmentation) . The insets show the exponential tail 
of the size distributions for shear fragmentation in the case of 
the moving convection cell flow and the flxed convection cell 
flow. 



average size classes match, it can be seen that the shape 
of the particle size distributions for both flows have be- 
come almost identical (Fig. 4) and show the characteris- 
tic exponential tail. Again, this indicates that the shape 
of the size distribution is mostly affected by fragmenta- 
tion and any flow specific differences in the distribution 
result mainly from differences in the shear distribution in 
the flow, which in turn change the fragmentation rates. 

This strong dependence of the steady state on the frag- 
mentation process can also clearly be seen by how ctoo 
changes with the coagulate strength 7 (see Fig. 5). We 
note that the shear forces in the convection flows and the 
sinusoidal shear flow have a different magnitude, as seen 
in Fig. I. Therefore three different values of the coagu- 
late strength parameter 7 need to be chosen to yield a size 
distribution that is not only determined by size-limiting 
fragmentation. For the convection flows the coagulate 
strength 7 needs to be approximately a factor 3 larger 
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FIG. 4: Histogram of the particle size distribution in steady 
state for the moving convection cell flow (7 = 50) and the 
smoothed out sinusoidal shear flow (7 = 8.5, (3 = 7r/20) for 
shear fragmentation. 



than for the sinusoidal shear flow. 

It is clear that Qqc increases with 7, because particles 
become more resistant to shear. The exact functional 
relationship is however not trivial. A flrst qualitative 
estimate of the shape of this aoo(7) curve can be derived 
by assuming that over one period of the flow the particles 
experience an 'average shear' 
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(7) 



where G(x, t) is the modulus of the local velocity gradi- 
ent, p(x, t) is the distribution of particles and D is the 
unit square domain. From Eq. (3) we then get for the 
average critical size at this velocity gradient 



= G-373 



(8) 



Particles that exceed this size will therefore typically 
break up during one period of the flow. Since particles 
break into two parts due to shear, the average size would 
then be aoo > acrit/S. The average shear G is, however, 
somewhat complicated to estimate. It would have to be 
calculated as a mean over the positions of all particles in 
the flow at a given time. Additionally, how much larger 
than the critical size particles get before they break up 
depends on the coagulation probabilities, and therefore 
also on the local concentrations of particles. The exact 
dependency of aoo(7) is therefore not easily calculated. 
What can be seen from Eq. (8) is however that the av- 
erage size is expected to scale with 7 as 



oc 7 



(9) 



This dependence is expected to hold for all values of 
7 and Q!oo, where shear fragmentation dominates. A fit 
with Eq. (9) for the different flows is shown in Fig. 5 
and for lower values of 7 the fits agree very well with the 
simulation results. It can be seen that for higher values 
of 7, when size-limiting fragmentation becomes impor- 
tant, the aoo(7) curves deviate from this estimate and 
converge towards the limiting value q;£™^ (see Fig. 5). 
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FIG. 5: Asymptotic average size class index Qoo as a function 
of the coagulate strength parameter 7. Squares: numerical 
results, solid line: fit by Ofoo — ci +£27^ (see Eq. (8)). (a) 
moving convection cell flow (ci = 0.521, C2 — 3.7 x 10^'') (b) 
fixed convection cell flow (ci = —0.256, C2 = 6.2 x 10^'') and 
(c) sinusoidal shear flow (ci = 0.945, C2 = 1.6 x 10"^). 



This result demonstrates how the steady state depends 
very strongly on the fragmentation process. However, 
the different proportionality constants for Eq. (9) still 
depend on the flow and can also depend on the spatial 
distribution of the particles, since different regions of the 
flow might exhibit different shear. 

Finally, we mention some further results from our 
model. First, for all flows the steady state in the case of 
shear fragmentation is not static, instead due to the pe- 
riodic time dependence of the flows the steady state also 
varies periodically over time. This is very clear for ex- 
ample for the fixed convection flows when looking at the 
average size class index {a{t)) = Ei=i "iA^a, 
where N^. denotes the number of particles in size class 
at (Fig. 6). 

Such a time dependence of the average particle size, 
and therefore of the whole particle size distribution, can 
have important physical consequences. For example the 
settling of coagulates in the ocean, which is an impor- 
tant part of the biological carbon pump in the ocean, 
is greatly affected by the size distribution of the coag- 




50 100 150 200 250 148 152 156 160 
t t 

FIG. 6: In the case of shear fragmentation the steady state 
is not static, instead it fluctuates with the period of the flow: 
(a) average size class vs. time for the flxed convection flow 
and (b) close-up of (a). 

ulates. In coastal areas, where the fluid may be peri- 
odically forced by the tides, such a time dependence of 
the distribution can greatly affect the settling rates and 
therefore the deposition of coagulates [47] . 

Second, since the range of possible parameter values for 
our model is very large, we also mention the robustness 
of our findings with respect to the different parameters: 



la) (b) (c) 
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FIG. 7: Normalized number of particles versus relative radius 
in steady state in the case of shear fragmentation for different 
values of the coagulate strength parameter 7. The normalized 
size distributions collapse onto a master curve that depends 
on the fragmentation mechanism and the shear forces in the 
flow for (a) moving convection cell flow (7 — 44 (circle), 50 
(square), 56 (triangle)), (b) fixed convection cell flow (7 = 
42.5 (circle), 45 (square), 50 (triangle)) and (c) sinusoidal 
shear flow (7 — 16 (circle), 17 (square), 18 (triangle), f3 = 
20/7r). 

(i) For each flow there is a certain range of the coag- 
ulate strength parameters 7 where the size distribution 
for shear fragmentation is "fully developed" . By this we 
mean that 7 is large enough so that a sufficiently large 
fraction of particles has left the smallest size class, but 
7 is small enough so that break-up due to size-limiting 
fragmentation does not play a significant role. In this 
intermediate 7 range, where the particle size distribution 
is fully developed, a scaling form 

^" = / (JL) (10) 

max(7V„) -^{{r)) ^ 
is found to hold, where (r) represents the average radius. 
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Note that the form of the size distribution is independent 
of 7. All distributions in this parameter range collapse 
then onto a single master curve. 

While this scaling form is independent of the param- 
eters of the coagulation and of the coagulate strength, 
the difference between the different flows remains. More 
specifically, the scaling form / changes when shear forces 
in the flow or the fragmentation mechanism, for example 
the distribution of fragments, is varied. 

(ii) When investigating cases with different total mass 
M, we find that for size-limiting fragmentation otao is 
largely independent of M . For shear fragmentation with 
M < 3x10^7711, aa^iM) increases approximately linearly 
with M, while for higher values a saturation of aoa{M) 
sets in, which is due to the fact that size-limiting frag- 
mentation dominates in this case. 

(iii) By considering other initial particle size distribu- 
tions than mentioned above, for example any single size 
class with a > 1 or a uniform distribution of sizes, and 
keeping the total mass M fixed, the asymptotic state is 
found to be independent of the chosen initial distribution 
for each flow and for both fragmentation rules. 

(iv) We also investigated the role of the number of size 
classes and found that in the chosen range of 7 values 
the size distributions for shear fragmentation are not in- 
fluenced by the number of size classes. 

(v) The effect of the number of new particles formed 
by fragmentation has been considered. For instance, the 
distributions of particles for ternary fragmentation are 
similar to the ones for binary splitting and only show a 
slight shift towards smaller size classes [49]. 

IV. SUMMARY 

We discussed the formation of a steady state size distri- 
bution in a coagulation fragmentation process of inertial 
particles transported by different flows. Our most im- 
portant finding is that fragmentation rather than coagu- 
lation is the dominating process for the steady state size 
distribution. For size-limiting fragmentation we found al- 



most no differences in the shape of the steady state size 
distribution for various flows. Even in flows with great 
differences in coagulation rates, e.g. due to differences in 
local particle concentrations, particle size distributions 
remained very similar. For the case of shear fragmenta- 
tion differences in the shape of the size distribution for 
the different flows appeared. It was found that these were 
due to differences in the spatial and temporal variation 
of the shear, which in turn affected the fragmentation. 

We have shown that an individual particle based mod- 
eling approach is able to reflect typical properties of coag- 
ulation and fragmentation processes of inertial particles. 
The appearance of a steady state is demonstrated. We 
outlined some of the differences in the convergence to the 
steady state and the particle size distribution that can re- 
sult from different types of fragmentation and flow. Al- 
together, our results suggest that coagulation dominates 
different time spans of the process than fragmentation. 
While coagulation is most important for the transients in 
the beginning, the steady state size distribution is mainly 
determined by fragmentation. As a consequence the spa- 
tial distribution of particles plays only a transient role. 
The underlying flow is important for the steady state 
in the case of shear fragmentation, as the spatial and 
temporal variation of the shear can greatly influence the 
fragmentation rates. 

The generalization to a fully three-dimensional system 
is straightforward. The relaxation towards the steady 
state would slow down, due to the decreased probabil- 
ity of collisions. However, our conclusions regarding the 
dependence of the steady state size distribution on the 
fragmentation remain valid. 
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